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A GLOBAL MAGNETIC ANOMALY MAP 
R. D. Regan*, J. C. Cain**, and W. M. Davis* 


ABSTRACT 

A subset of POGO satellite magnetometer data has been 
formed that is suitable for analysis of crustal magnetic anon- 
malies. Using a thirteenth order field model, fit to these 
data, magnetic residuals have been calculated over the 
world to latitude limits of ±50°. These residuals averaged 
over one degree latitude -longitude blocks represent a de- 
tailed global magnetic anomaly map derived solely from 
satellite data. Preliminary analysis of the map indicates 
that the anomalies are real and of geological origin. 
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A GLOBAL MAGNETIC ANOMALY MAP 


INTRODUCTION 

The U. S. Geological Survey and Goddard Space Flight Center (NASA) are jointly 
involved in an analysis of presently available magnetic data from the Cosmos-49 
(Benkova et al. , 1971) and OGO-2, 4, and 6 (POGO) (Cain and Langel, 1971) 
satellites. The major objective of this analysis is the identification of anomalies 
having geological significance. 

Previous reports have established the gross correlation between aeromagnetic 
surveys and satellite measurements over broad areas and that magnetic anoma- 
lies associated with inhomogenie ties in the lithosphere should be detectable at 
satellite altitudes (Zietz et al. , 1970; Andreasen et al, , 1970). The first con- 
firmation of the existence of such an anomaly in the satellite data was the identi- 
fication of the Bangui anomaly (Regan et al. , 1973a). This anomaly, occuring 
over central Africa, was discovered in a detailed analysis of the POGO data and 
later identified on Project MAGNET (Stockard, 1971) data. 

The data reduction procedures employed to define the Bangui anomaly have now 
been applied to the entire set of POGO magnetic measurements. It was thus 
possible to produce a detailed global magnetic anomaly map. The resolution of 
this map permits the identification of anomalies with wavelengths less than 14°. 
This represents a substantial increase in resolution over previous analyses 
(Benkova et al. , 1973; Regan et al. , 1973b). The anomalies noted on these 
earlier maps termed "intermediate size" anomalies, were of half wavelengths 
of approximately 19° and not clearly of lithospheric origin. 


DATA SELECTION AND ANOMALY CALCULATION 

At satellite altitudes, the anomaly signal is comparable in magnitude to time 
variations from external sources, therefore the satellite data must be obtained 
when the field is undisturbed. Also, the rapid decay of the signal with altitude 
requires that measurements be obtained at as low an altitude as possible. To 
obtain such data, the POGO observations, made during a period of variable 
magnetic activity and over an altitude range of 400 to 1500 km, had to be screened 
for low altitude and minimal magnetic activity. 

A subset of the POGO data was constructed by screening the observations on 
altitude, local time, and Planetary Magnetic Activity Index (Kp) (Sugiura and 
Heppner, 1968, p. 72). Only data collected when Kp was less than, or equal to 
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2 + were utilized. These data were corrected for disturbance effects by the 
method described by Davis and Cain (1973). Because of the limitations of this 
disturbance correction, the data were further restricted to geographic latitudes 
less than, or equal to ±50°. Of this set, data collected at altitudes greater 
than 700 km and at local times between 0900 and 1500 were rejected. The local 
time selection was made to avoid magnetic effects due to the equatorial electro- 
jet and ionospheric currents (Cain and Sweeney, 1973). This screening of the 
data provided a subset of 393,452 observations with a 7 second or 50km spacing 
along a given orbit and good geographic distribution (typically 15 observations 
per 1° latitude -longitude block). 

The method utilized to determine an anomaly is a regional-residual field sepa- 
ration where the regional field is defined by a four-dimensional (latitude, longi- 
tude, altitude, time) function (Cain et al. , 1967). This function is computed by 
least squares fitting of a spherical harmonic series to the observed data. The 
series is truncated at a harmonic order where the root mean square of the fit 
does not decrease significantly by computing higher order coefficients or at some 
prior arbitrary point necessitated by computer capacity. The anomalies or 
residuals, termed AF, are defined as the difference between the measured and 
computed values at each observation. 

The validity of field residuals at satellite altitudes has been determined by com- 
paring residual values derived from two completely independent sets of obser- 
vations and field references (Regan et al. , 1973b). 


WORLD MAP 

A thirteenth order spherical harmonic series (or field model) was fit to the sub- 
set of POGO data. Magnetic field residuals were calculated and averaged over 
one degree latitude -longitude blocks. The resulting global magnetic anomaly 
map, contoured at 2 gammas, exhibited some north-south striping that was a 
reflection of the direction of the POGO orbits. A preliminary attempt to reduce 
some of the noise in the map was accomplished by deleting the ±2 gamma contours 
that were not associated with more intense anomalies; that is, only anomalies 
with absolute magnitude greater than 2 gammas were retained. This filtered 
map is shown in Figure 1. This map represents one degree averaged magnetic 
residuals at a mean elevation of 520 km. The residuals are not at a common 
elevation, but on a topographic surface defined by the averaged elevations of 
the data in each one degree block. However, the anomalies are not considered 
to be significantly distorted by this variation in altitude. Only minor changes in 
anomaly structure were noted in the reduction to a common plane of a 20 by 20 
degree area of this map over the Bangui anomaly (Regan et al. , 1973a). 
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The maximum and minimum anomaly values are +10 and -12 gammas respec- 
tively. Colors were utilized to simplify the complex structure of the map. 

Color changes occur at the contour interval of two gammas. The red end of the 
spectrum is associated with positive anomalies and the blue end with negative 
anomalies. Specifically, the colors used in this figure are: purple, -4 gammas 
and below; dark blue, -4 to -2 gammas; light blue, -2 to 0 gammas; green, 0 to 
+2 gammas; yellow, 2 to 4 gammas, orange, 4 to 6 gammas; pink, 6 to 8 gam- 
mas, and red, 8 gammas and above. 

Several features are immediately apparent in the global magnetic anomaly map 
{Fig. 1). Most striking is the apparent lack of anomalies in the Pacific Ocean 
with the exception of a pair of anomalies (+6 gammas and -6 gammas amplitude) 
in the area where the Emperor Seamount chain and the Hawaiian Island chain 
intersect. Also evident is the Bangui anomaly, a 12 gamma low, over central 
Africa. 

The Bangui anomaly is a confirmed magnetic anomaly of possible lithospheric 
origin. Because of this it is presumed that the other magnetic anomalies shown 
are also real and of internal origin. However, it is possible that a few of the 
anomalies could be artifacts of spurious orbits, insufficient time corrections, 
and other data reduction procedures. Thus, the validity of all the anomalies is 
now being investigated by several approaches. 

First the anomaly structure as revealed on individual satellite passes is being 
examined. The anomalies shown on individual passes are not affected by the 
averaging technique used to construct the map. Also, persistence of an anomaly 
in satellite observations made at various local times and at various elevations 
would indicate that the anomaly was real and not the result of any external field 
effects. The magnetic field residuals on several POGO satellite passes over 
India are shown in Figure 2. The magnetic high over the southern part of India 
and the low over the Himalaya Mountains are shown to persist on all satellite 
passes. The anomalies are independent of local time and show a dependence on 
altitude that is consistent with an internal origin. 

Another method of investigating the validity of these anomalies is by comparing 
them to Project MAGNET (aeromagnetic) data. The available Project MAGNET 
data over several anomalous areas of the satellite map have been examined and 
in each case there is evidence of an anomaly. The signals associated with the 
Bangui magnetic anomaly have been compared by continuing the Project MAGNET 
data to satellite altitude, as shown in Figure 3. This represents a segment of 
Project MAGNET Line T217 over the Bangui magnetic anomaly. Assuming a 
two-dimensional source the profile was mathematically continued to satellite 
altitude and compared to a satellite profile. The comparison demonstrates that 
the same anomaly is being measured by both satellite and aircraft magnetom- 
eters. The magnitude of the anomaly at aircraft elevation (-780 gamma) 
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compared to that at satellite elevation (-16 gammas) is also consistent with a 
magnetic anomaly of internal origin. 

An alternative approach is to compare the POGO magnetic anomalies with a 
similar reduction of the Cosmos-49 data. Such a comparison at longer wave- 
lengths was made by Regan et al, , (1973b) and demonstrated a correlation be- 
tween the two data sets. However, it was noted then that the Cosmos data are 
much less numerous than the POGO observations. Also, the limited orbital 
information in the Cosmos data is a source of error in the reduction of the 
magnetic measurements. Hence it was not thought that any meaningful results 
are avilable with a one degree averaging of the Cosmos data. Nevertheless a 
thirteenth order field model was fit to the Cosmos-49 data. As expected, the 
scatter in this fit was much higher (38 gammas rms) than that of the POGO 
model (7 gammas rms), and the data set was not dense enough for averaging 
over one degree latitude -longitude blocks. However, a map of the Cosmos 
residuals using this model averaged over five degree latitude -longitude blocks 
was constructed. Although the five degree averaged map is not as detailed as the 
POGO one degree averaged map (Fig. 1) the general location of the anomalies 
on both maps is consistent. An example is shown in Figure 4. This shows the 
Cosmos-49 data over the African continent. The POGO data, averaged on one 
degree latitude -longitude blocks, over the same area is shown in Figure 5. In 
general there is a good correlation between the two maps, particularly over the 
southern half of the continent. The distortions in the Cosmos map are primarily 
a result of the sparse nature and nonuniform distribution of the data. The num- 
ber of Cosmos observations per five degree block in this area ranges from 1 to 
23. Whereas the number of POGO observations ranges from 300 to 450 per five 
degree block. There is also a distortion of the anomalies along the flight path 
of the Cosmos satellite (inclination ~50 degrees), similar to the north-south 
striping seen on the POGO map (Fig. 1). 

Origin of the Anomalies 


It is suspected that the anomalies shown in this map are primarily of geologic 
origin with a source in the lithospheric region of the earth. 

The occurrence of these anomalies on all satellite passes, independent of local 
time and their decay with altitude, imply a definite internal origin for the source 
of these anomalies. Due to the use of a thirteenth order field model, the anomaly 
half wavelengths are less than 14° (~180°/13) (Kaula, 1967) and this shortwave- 
length end of the geomagnetic spectrum has been attributed to the crustal com- 
ponent of the geomagnetic field (Alldredge et al. , 1963). Also by considering the 
rate of decay of anomalies described by the fourteenth harmonic, an estimate of 
the magnetic intensity of the source relative to the satellite anomalies can be 
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calculated. A source in the lower mantle or at the mantle-core boundary would 
require an amplitude several thousand times as intense, whereas a source in 
the upper mantle or crust would need to be only an order of magntiude more 
intense. 

A second approach to the identification of the satellite magnetic anomalies as 
lithospheric is to examine their correlation with known tectonic, geologic, or 
geophysical data. This is now being done in detail for the African continent. 

An example is shown in Figure 5. Here a tectonic map of Africa {from Wyllie, 
1971), is shown with the POGO satellite magnetic data. The color scheme for 
the tectonic map is the same as that shown in Figure 6. Note that the color 
scheme for the magnetic map differs slightly from that of Figure 1. In this map 
the contour interval is 2 gammas and purple represents -6 gammas and below; 
dark blue, -6 to -4 gammas; light blue, -4 to -2 gammas; green -2 to -0 gam- 
mas; yellow, 0 to +2 gammas; orange, 2 to 4 gammas; pink, 4 to 6 gammas; 
and red, 6 gammas and above. 

In considering the magnetic anomalies in this figure and Figure 6, it must be 
noted that the inclination of the main field is varying as a function of latitude. 
Thus, anomalies, due to induction, having similar origins would have different 
magnetic signatures at various geomagnetic latitudes. 

In the center of the magnetic anomaly map is the Bangui anomaly that occurs 
over the tectonic uplift zone between the Chad and Congo Basins. South of the 
Bangui anomaly there is a general correlation of magnetic highs with continental 
shield areas and magnetic lows with continental platforms. The low at the 
southern tip of the continent is coincident with the Karroo Basin. The correla- 
tion appears reversed on the other side of the magnetic equator. The major 
magnetic high in the northwest part of Africa is situated over the El-Juf or 
Taudine Basin and the magnetic low east of this occurs over the uplift area be- 
tween this basin and the Chad Basin. 

For the general correlation on a global scale the filtered map (Fig. 1) is shown 
overlain on a tectonic map of the world (from Wyllie, 1971) in Figure 6. Only 
the contour lines from the magnetic map are shown. Figure 1 will have to be 
examined to determine the sense of the anomalies. Some correlations between 
the tectonic and magnetic maps are obvious. The continental shields in the 
southern hemisphere, and India in the north, all have magnetic highs associated 
with them. There is a considerable amount of magnetic structure in the vicinity 
of oceanic ridges contrasted with the central Pacific region that is devoid of 
anomalies. 
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CONCLUSION 


Analysis of the POGO satellite data has resulted in a detailed global magnetic 
anomaly map. Verification of several distinct anomalies was obtained by ex- 
amining individual satellite profiles and Project MAGNET data. The persistence 
of the anomalies on all satellite passes, independent of local time, demonstrate 
that they are real and not the effect of any magnetospheric or ionospheric dis- 
turbance. The rate of decay of the anomalies with altitude as observed in the 
aircraft and satellite data imply an internal origin for these anomalies. The 
limited analysis of the Cosmos data reveals a similar global distribution of 
anomalies which are more intense. This is consistent with an internal origin 
because the Cosmos-49 satellite was at a lower mean elevation ( ~375km) than 
the POGO satellites. Although much more work is needed to fully interpret the 
map and determine the causes of these anomalies, the initial comparisons with 
tectonic maps are quite encouraging. 

These results further demonstrate the utility of a satellite magnetometer as a 
geological/ geophysical tool. A proposed global survey at low altitude using a 
vector magnetometer would provide the increased resolution necessary for more 
details d geological dose ription. 
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Figure 1. Global Magnetic- Anomaly Map; Contour Intervals, 2 gammas. 
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Figure 2. Magnetic Field Residuals from Individual POGO Satellites. Passes 
over India; Note the magnetic high over Southern part of India and the low over 
Himalaya Mountains. 
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Figure 3. 


A Comparison of Bangui Magnetic- Anomaly at OGO 
Satellite and Aircraft Altitudes. 


12 






</> 

< 

2 

2 

< 

O 

U. 

<3 



Figure 4. Five Degree Averaged Cosmos-49 Magnetic Anomaly Map over 

Africa 
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Figure 5. Tectonic and POGO Satellite Magnetic- Anomaly Maps Over 
Africa (tectonic mapfrom Wyllie, 1971); Contour Intervals, 2 gammas. 
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Figure 6. Global Tectonic Map Combined with the Contours of the Satellite 
Magnetic-Anomaly Map (tectonic map from Wyllie, 1971); Contour Inter- 
vals, 2 gammas. 
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